home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slarf.z / slarf
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAARRRRFFFF((((3333FFFF))))                                                            SSSSLLLLAAAARRRRFFFF((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLARF - applie a real elementary reflector H to a real m by n matrix C,
  10.      from either the left or the right
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SLARF( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
  14.  
  15.          CHARACTER     SIDE
  16.  
  17.          INTEGER       INCV, LDC, M, N
  18.  
  19.          REAL          TAU
  20.  
  21.          REAL          C( LDC, * ), V( * ), WORK( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      SLARF applies a real elementary reflector H to a real m by n matrix C,
  25.      from either the left or the right. H is represented in the form
  26.  
  27.            H = I - tau * v * v'
  28.  
  29.      where tau is a real scalar and v is a real vector.
  30.  
  31.      If tau = 0, then H is taken to be the unit matrix.
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      SIDE    (input) CHARACTER*1
  36.              = 'L': form  H * C
  37.              = 'R': form  C * H
  38.  
  39.      M       (input) INTEGER
  40.              The number of rows of the matrix C.
  41.  
  42.      N       (input) INTEGER
  43.              The number of columns of the matrix C.
  44.  
  45.      V       (input) REAL array, dimension
  46.              (1 + (M-1)*abs(INCV)) if SIDE = 'L' or (1 + (N-1)*abs(INCV)) if
  47.              SIDE = 'R' The vector v in the representation of H. V is not used
  48.              if TAU = 0.
  49.  
  50.      INCV    (input) INTEGER
  51.              The increment between elements of v. INCV <> 0.
  52.  
  53.      TAU     (input) REAL
  54.              The value tau in the representation of H.
  55.  
  56.      C       (input/output) REAL array, dimension (LDC,N)
  57.              On entry, the m by n matrix C.  On exit, C is overwritten by the
  58.              matrix H * C if SIDE = 'L', or C * H if SIDE = 'R'.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAARRRRFFFF((((3333FFFF))))                                                            SSSSLLLLAAAARRRRFFFF((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDC     (input) INTEGER
  75.              The leading dimension of the array C. LDC >= max(1,M).
  76.  
  77.      WORK    (workspace) REAL array, dimension
  78.              (N) if SIDE = 'L' or (M) if SIDE = 'R'
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.